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ABSTRACT 

We have measured the Sunyaev-Zcldovich effect towards the Coma cluster (A1656) with the MITO 
experiment, a 2.6-m telescope equipped with a 4-channel 17 arcminute (FWHM) photometer. Measure- 
ments at frequency bands 143±15, 214±15, 272±16 and 353±13 GHz, were made during 120 drift scans 
of Coma. We describe the observations and data analysis that involved extraction of the SZ signal by 
employing a spatial and spectral de-correlation scheme to remove a dominant atmospheric component. 
The deduced values of the thermal SZ effect in the first three bands are ATo=— 179±38,— 33±81,170±35 
/jK in the cluster center. The corresponding optical depth, r=(4.1±0.9)xl0~ 3 , is consistent (within 
errors) with both the value from a previous low frequency SZ measurement, and the value predicted 
from the X-ray deduced gas parameters. 

Subject headings: cosmology: cosmic microwave background - observations - galaxies: clusters: 
individual (A1656) 



1. INTRODUCTION 

Compton scattering of the cosmic microwave back- 
ground (CMB) radiation by electrons in the hot gas in 
clusters of galaxies - the Sunyaev-Zeldovich (SZ) effect - 
has long been recognized as a uniquely important feature, 
rich in cosmological and astrophysical information (Zel- 
dovich & Sunyaev 1969; Sunyaev & Zeldovich 1970; Boyn- 
ton & Melchiorri 1978; Gunn 1978; Silk & White 1978; 
Cavaliere et al. 1979). 

Interest in the effect heightened when high quality 
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images of the effect were obtained with interferometric 
arrays (Carlstrom et al. 2001; for general reviews, see 
Rephaeli 1995a, Birkinshaw 1999). 

Accurate determinations of the Hubble constant (Ho), 
cluster mass profiles, and cluster peculiar velocities from 
SZ and X-ray measurements require precise description of 
the intracluster (IC) gas temperature and density profiles, 
and improved control of systematic errors. These can be 
more optimally achieved by multi-frequency SZ measure- 
ments, and high-quality spectral and spatial X-ray mea- 
surements of nearby clusters. 

Therefore, it is essential to develop sensitive SZ experi- 
ments operating at telescopes with effective beam sizes of 
a few arcminutes, equipped with bolometer arrays oper- 
ating at frequencies on both sides of the Planckian peak. 
One such ground-based system is MITO (Millimetre and 
Infrared Testagrigia Observatory; De Petris et al. 1996). 

Coma is the richest nearby cluster whose IC gas prop- 
erties have been determined from extensive X-ray obser- 
vations. Attempts to measure the S-Z effect in Coma at 
low frequencies were reported by Pariskii (1973), Rudnick 
(1978), Lake & Partridge (1980), Birkinshaw et al. (1981), 
and Herbig et al. (1995, revised in Mason et al. (2001)) 
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while at high frequencies by Silverbcrg et al. (1997). 

2. OBSERVATIONS 

2.1. MITO 

MITO is a 2.6-m alt-azimuthal telescope located at an 
altitude of 3480 m in the Italian Alps. The telescope de- 
sign, optical configuration, and electro-mechanical modu- 
lation system are described in Dc Pctris et al. 1989, Ger- 
vasi et al. 1998, Mainella et al. 1996. The 4-channel 
photometer (Orlando et al. 2002) consists of 4 NTD Ge 
composite bolometers cooled down to 290 mK. The ther- 
modynamic NET is around 1 mK s 1 / 2 . At the telescope 
focal plane the photometer has a 17' FWHM field of view 
(FOV). 

Observations were made during February and March of 
2000 & 2001, when data from some 120 drift scans of Coma 
were collected. Each scan consists of 10 minutes of mea- 
surements collected every second. A 41' beamthrow was 
selected for a 3-field square-wave like spatial modulation. 
The position of Coma in the sky covered altitude values 
from 50 up to 73 degrees above the horizon. The drift scan 
systematics are better controlled than in source tracking 
mode, mainly stemming from weaker microphonics and 
lower changes in sidelobe spillover effects. At a data sam- 
pling rate of 1 Hz, pointing information, signals, modula- 
tion system and photometer performance were recorded. 
A pointing accuracy of ~ 1' was attained by frequent mea- 
surements of bright stars close to the source by a CCD 
camera. We measured beam shape, responsivity and at- 
mospheric transmission for each channel by means of plan- 
etary calibrators - Jupiter, Saturn and the Moon - before 
each Coma observation session. In order to correct for 
atmospheric absorption we measured the channels zenith 
transmission during each night. 

2.2. Results from X-ray Observations 

Coma has been extensively observed, mostly in the 1- 
10 kcV range. The most recent high spectral and spa- 
tial resolution measurements with the XMM/EPIC ex- 
periment are best fit by fcT=8.2±0.1 keV (at 90% confi- 
dence level) in the central 10' region (Arnaud et al. 2001). 
Analysis of the emission from a larger region of the clus- 
ter yields fcT=8.2±0.4 keV, which we adopt here. This 
value of the temperature is in very good agreement with 
the value deduced from Ginga observations. Since quan- 
titative results for the gas density profile from neither 
XMM nor Chandra are currently available, we use the 
ROSAT (Briel, Henry, & Bohringer 1992) deduced values 
n e0 ~ (2.89 ± 0.04) x 10" 3 cm" 3 , r c ~ 10.5' ± 0.6', and 
(3 ~ 0.75, for the central electron density, core radius, and 
index in the expression for the commonly used f3 density 
profile, n e (r) = n e0 (l + r 2 /r 2 )~ 3/3 / 2 . 

3. DATA ANALYSIS 

Clearly, there is no unique approach to the analysis of 
a noisy dataset which includes a dominant, fluctuating at- 
mospheric emission, and other confusing signals. When 
the spatial distribution of the confusing source is unknown, 
such as that of CMB anisotropy, a sufficient number of 
channels is selected so as to allow removal of anticipated 
foregrounds but the corresponding system of equations can 
often be numerically unstable. 



Our task is made feasible largely by the known posi- 
tion of the source, and its approximate angular extent. 
Therefore, to first approximation, we can fit the data 
with a signal of a known shape but unknown peak am- 
plitude. The errors in measuring this amplitude approach 
the intrinsic detector noise if we can successfully remove 
the contributions of atmospheric fluctuations and CMB 
anisotropy whose spectral shape is known. In addition to 
the SZ effect, these signals were modeled separately from 
all others which were collectively treated as unidentified 
noise. While three channels are in principle sufficient to 
achieve this goal, we added a fourth channel centered on 
the crossover frequency (where the thermal SZ effect van- 
ishes) , in order to check the separation between CMB and 
atmospheric emission, and to possibly obtain a rough esti- 
mate on the value of the kinematic SZ effect. 

We have computed the expected spectral ratios of the 
SZ signals in channels 2, 3, and 4 with respect to the 
first channel. The calculation is relativistically exact and 
is based on the treatment of Rcphacli (1995b; see also 
Rephaeli & Yankovitch 1997), taking the above quoted 
range of the gas temperature. At each frequency we cal- 
culated the intensity change due to the thermal SZ effect, 
AI SZ , and convolved it with the spectral response of the 
photometer, £j, and the atmospheric transmittance e a tm, 
which was determined from a model (Liebe 1985). The 
expected ratio of the i-th to the first channel is thus, 

b _ AiliRi / °° AI sz (v)ei(is)e atm (is)dv 
Af2ii?i f °° M sz {v)e- ^ {v)e a tm{v)dv , 

where Ri is the responsivity in the i-th channel as mea- 
sured with sky calibrators and A£li is the throughput, 
which slightly changed from channel to channel due to 
the different optical paths through the telescope and pho- 
tometer. Uncertainties in the evaluation of the bi are due 
to the relatively small error in the value of kT and the 
fluctuating values of atmospheric transmittance (0.5-2.0 
mm of precipitable water vapour). Similarly, we have cal- 
culated the ratios, Ci, of the intensity change due to the 
primary CMB anisotropy in the four bands. In Table 1 all 
the cited ratios are listed with Ri and AVLi. 

The recorded data in the four channels (i=l,..,4) for each 
scan (j = l,..,120) can be written as follows: 

AS? = a?Al& m (i) + kAV£(M) 

+ Ci AVi mb (a,5) + AVi ffset .(a,5;t) 

+ W s j pike .(t) + AVi oise4 (t)- (2) 

In this equation, biAVj z is the SZ signal which is fit 
with wlAV s 3 z _ sm , where AV s 3 z _ sm is the simulated SZ 

signal, normalized to 1 at the peak; u>l are the esti- 
mated peak values of SZ signals (in the absence of noise 
= b^ 1 (wi),i = 2,3,4) and al are the atmospheric ra- 
tios with respect to channel 1, to be determined for each 
scan. The al values could change from scan to scan due 
to the atmospheric long time variations while they are as- 
sumed to be constant within a single scan. Since atmo- 
spheric fluctuations dominate on each scan, a rough esti- 
mate of the al can be done as the ratios of the channels 
standard deviations. For this purpose we used only the 
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data outside the expected position of Coma. In Table 1 
we have also quoted the mean values of a, averaged over 
all the drift scans in each observational campaign. Detec- 
tor and other sources of noise are collectively included in 
the term AV^ „• ,. Estimated errors in the evaluation of 

w\ due to uncertainties in hi constitute a small fraction of 
the overall error. 

The data include three main sources of contamination: 

a. Cosmic ray spikes (SV^ pike we observed roughly 
one spike in every 10 scans. A special de-glitching algo- 
rithm has been applied to correct these data. However, 
their total removal from the analysis did not change our 
results apart from a slight increase of the noise. After re- 
moval of the spikes the data were averaged over 15 seconds 
so that each drift scan consists of 37 bins roughly separated 
in sky by 3.75 arcminutes, about 1/4 of the FOV. 

b. An instrumental offset (AV^j: set i (a 1 d; t)) ranging 
from 0.4 /iV (channel 1) up to 5/iV (channel 4) was mea- 
sured. This is partly due to the residual side lobes of the 
telescope, and partly to the thermal imbalance of the pri- 
mary mirror. The telescope is stopped during each drift 
scan. In this way the offset is obviously independent of po- 
sition, so it is possible to remove it from each scan through 
a linear fit of the data; a quadratic fit did not change 
the rms value of the residual fluctuations. Note that the 
removal of the offset through a linear fit may introduce 
an additional effect, due to the presence of large atmo- 
spheric fluctuations. The average of these along a finite 
drift scan will introduce an additional offset, AV^j_ atm t 

= -gj Y^kLi a \^Vitm k (*)■ This unknown offset will in turn 
introduce an error in the estimated value of the SZ signal; 
we discuss below how this offset was removed. 

c. Large atmospheric fluctuations; we based our approach 
to minimize their effect on the measured data in channel 
4, based on the fact that a 4 is ~ 13. The efficiency of this 
procedure depends on the amplitude of the atmospheric 
fluctuations in each drift scan, as well as on the degree 
of correlation among the various channels with channel 
4. A detailed discussion of the atmospheric noise will be 
given elsewhere; here we briefly discuss only the essen- 
tial points. In channel 2 the residuals of the signal after 
subtracting a^/aP^AS^ and minimizing the difference con- 
sist of uncorrected atmospheric fluctuations, signal due to 
CMB anisotropy, and only a very small SZ signal. If we 
fit the residuals with a simulated SZ signal of amplitude 
tiP 2 i we S e t (^2) — —0.25 ± 0.32 nV. This result indicates 
that the SZ signal in channel 2 is consistent with zero, as 
expected, and the large dispersion - with respect to the 
final error due to the intrinsic detector noise of about 0.05 
nV - around the mean suggests that there is a residual 
atmospheric contamination Aw2 = a 2 Awi. A significant 
decrease in the dispersion is obtained if we remove the val- 
ues of W2 far from the mean. These values are possibly due 
to the presence of a large atmospheric disturbance along 
the line of sight (los) to Coma. For instance, if we remove 
values which are more than 4a (standard deviation) from 
the mean, we obtain (u^) — —0.04 ± 0.09 nV. A similar 
analysis has been done for channels 1 and 3, for which 
we expect SZ signals at the predicted ratio 63. Residual 
atmospheric fluctuations - as well as the spurious offset 
discussed previously - introduce an additional error. The 



corresponding small corrections Awi and Aw^ — a^Awi, 
which have to be added to the measured values of the SZ 
effect in order to get the correct value, must satisfy the 
relation: w\ + Aw\ = (u> 3 + a^Awij/b^. The dispersion 
of (u>i + Au>i) is minimized by changing 63 within ±1.5 
63. The dispersion decreases by a factor of four when 63 
coincides with 63 within the errors. 

After correcting for binning and atmospheric transmit- 
tance we obtain (w*) ~ -0.171 ± 0.036, (w*) ~ -0.037 ± 
0.090, and (w*) ~ 0.174 ± 0.036 nV. In Figure 1, we show 
plots of AV SZ , AV cm b 1 and AV atm as derived by solving 
eq.(2). The data have been filtered with an adaptive filter 
which removes spatial patterns that are not present in the 
expected SZ profile. The filter does not alter the SZ signal, 
but strongly reduces CMB and atmospheric fluctuations. 




Fig. 1. SZ, CMB, and atmospheric signals extracted from 
equation (2), corresponding to the 143 GHz band, averaged 
over the drift scans and binned in 9 independent sky posi- 
tions. Atmospheric fluctuations have an amplitude compara- 
ble to that of the SZ effect, even after 20 hours of integration. 
CMB anisotropics are averaged over various parallactic angles 
corresponding to the different drift scans. 

The w value has been converted to thermodynamic 
temperature for the i-th channel by calibration based on 
measurements of Jupiter (Encrenaz & Moreno 2002) . In- 
tegrating its spectral temperature distribution over the 
four channels we determined the optical responsivities. 
The overall calibration error is about 10%, dominated by 
the uncertainty in the brightness temperature of Jupiter 
(~ 5%), and statistical errors in the estimate of the 
beamsize. The final values of the deduced temperature 
changes are AT ~ -73.5±15. 5,-13. 7±33.3,69.6±14.4 ^K, 
in the first three channels, respectively. Based on our 
calculated value for the geometrical form factor due to 
beam dilution and modulation strategy, 77 = 0.41±0.02, we 
compute the corresponding values AT ~ -179. 3± 37.8,- 
33.4±81.2,169.8±35.1 for AT in the cluster center. 
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4. DISCUSSION 

MITO measurements of the Coma are the first successful 
detection of the SZ effect in this cluster at high frequen- 
cies. Our results imply a mean Thomson optical depth 
r = (4.1±0.9) x 10~ 3 for a los through the center of Coma. 
Herbig et al. (1995) have already detected the effect at 
the much lower frequency of 32 GHz. Their reported cen- 
tral value, AT ~— 505±92 /iK, and corresponding optical 
depth, r = (5.6± 1.1) x 10~ 3 , were recently updated (Ma- 
son et al, 2001) with new calibration data, yielding an ob- 
served decrement AT—— 302±54 fiK. This value, corrected 
for the dilution over their 7.3' FOV and the switching am- 
plitude of 22.2', corresponds to AT =-520±93 fiK and 
t = (4.0 ± 0.6) x 10~ 3 , and is clearly in good agreement 
with X-ray deduced gas parameters. Clearly, these two 
results for To arc quite consistent, given the appreciable 
(la) errors. Even though our somewhat lower value is 
not significantly different from that of Herbig et al., it is 
interesting to note that the CMB anisotropy signal that 
we seem to have identified along a los to Coma (heav- 
ily exploiting our multi-frequency capability) could have 
possibly contributed to the SZ signal at the low (single) 
frequency OVRO measurement (where the two signals can- 
not be spectrally separated). Irrespective of whether this 
is indeed the case, the great advantage of multi-frequency 
work has been clearly demonstrated in the analysis of our 
measurements. 

In Figure 2 we show the first SZ spectrum of Coma by 
combining OVRO with our measurements; a fit to the pre- 
dicted spectrum yields t ~ (4.9 ± 0.7) x 10~ 3 . 

In addition to measurement errors, other errors include 
the use of a simple one dimensional, isothermal gas density 
with a (3 profile. Based on the X-ray morphology of Coma, 
the relatively small degree of cllipticity of X-ray morpholo- 
gies of other rich clusters, and the recent verification of the 
gas isothermality in Coma from XMM measurements (Ar- 
naud et al. 2001), we estimate that the additional error 
introduced by taking the gas to be spherically symmetric 
and isothermal is < 30% of the overall quoted error. Note 
also that the additional Comptonization by non-thermal 
electrons, (whose existence in Coma is deduced from mea- 
surements of extended radio emission) is negligible (Shi- 
mon & Rephacli 2002, Colafrancesco, Marchegiani & Pal- 
ladino 2001) if the energy density in these electrons is at 
the level deduced from measurements of high energy (> 20 



keV) X-ray emission with the RXTE (Rephaeli, Gruber & 
Blanco 1999) and BeppoSAX (Fusco-Femiano et al. 1999) 
satellites. This result is valid also for other viable energetic 
electron models (Shimon & Rephacli 2002). 
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FIG. 2. SZ spectrum of the Coma cluster. The solid 
line shows the best fit spectrum (taking isothermal gas with 
kT = 8.2 keV) to the combined MITO and OVRO (Mason et 
al. 2001) measurements, corresponding to r ~ (4.9±0.7) x 10~ 3 . 

The work reported here is only the first stage of an ex- 
tensive project to upgrade MITO to a bolometer array con- 
sisting of 36 elements operating at four frequency bands 
and with beamsizes down to 4'. Work on the upgraded 
system is underway, with first observations scheduled for 
early 2003. We plan to measure the effect in a large sam- 
ple of nearby clusters, with the goal of determining H and 
cluster masses. Based on the results from our long cam- 
paign to measure the effect in Coma, we are reasonably 
confident that this project is quite feasible. 

Contributions to the success of the MITO program were 
made by many scientists; in particular, we are grateful to: 
S. Masi, who designed and built the cryogenic He-3 sys- 
tem; P. Lubin, who kindly supplied a low noise cryogenic 
FET; T. Encrcnaz, for her FIR spectra of planets, and to 
P. de Bernardis for various useful suggestions. This work 
has been supported by CNR (Testa Grigia Laboratory is 
a facility of Sezione IFSI in Turin) COFIN-MIUR 1998, & 
2000, by ASI contract BAR, and by a NATO Grant. 
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